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Abstract 



We introduce semiclassical methods into the study of the volume spectrum in loop gravity. 
The classical system behind a 4-valent spinnetwork node is a Euclidean tetrahedron. We investi- 
gate the tetrahedral volume dynamics on phase space and apply Bohr-Sommerfeld quantization 
to find the volume spectrum. The analysis shows a remarkable quantitative agreement with the 
volume spectrum computed in loop gravity. Moreover, it provides new geometrical insights into 
the degeneracy of this spectrum and the maximum and minimum eigenvalues of the volume on 
intertwiner space. 

1 Introduction 

Knowing the classical geometry hiding behind a quantum system provides many insights into its 
behavior pQ. In this paper we introduce semiclassical methods to investigate the main feature of 
loop gravity: the prediction of a quantum discreteness of the volume of space [21 [3] . We focus on an 
elementary 4-valent grain of space and derive the semiclassical spectrum of the volume exploiting 
knowledge of the classical system associated to it — a Euclidean tetrahedron. This analysis devel- 
ops the Bohr-Sommerfeld quantization of space introduced in [2] and provides a new approximation 
scheme that can be tested against the numerical results of [21 El El El [101 E] • Moreover, the semi- 
classical geometries that arise give new insights into the scaling of the spectrum in a variety of limits. 

In loop gravity |12} [T3| HH] quantum states are built on a spin colored graph. The nodes of this 
graph represent three dimensional grains of space and the links of the graph encode the adjacency 
of these regions. Traditionally [151 ES] the quantum volume of a region of space R is obtained by 
regularizing and quantizing the classical expression 



where h is the metric of space. The total volume is obtained by summing the contributions from each 
node of the spin colored graph contained in the region R. Due to the necessity of a regularization 
scheme there are different proposals for the volume operator at a node. The operator originally 
proposed by Rovelli and Smolin [2] is 




(1) 
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where E r are electric-flux operators associated to each link of the graph impinging on the node, 
a > is a multiplicative constant, and the sum over r < s < t spans triples of links at the node. A 
second operator introduced by Ashtekar and Lewandowski [3] is 



Val = a \ V e(v r ,v s ,vt) E r ■ (E s x E t ) , (3) 




where v r are the tangents to the links of the spin network graph T embedded in the space 3- 
manifold. The e are signs, e(v r ,v s ,v t ) = ±1 or corresponding to whether the triple of tangents 
is right handed, left handed or planar respectively. When Val is viewed as an operator on the 
intertwiner Hilbert space associated to the node, the tangents v r and the signs e(y r , v s ,v t ) have to be 
understood as external fixed data. This residual dependence on the embedding of the spin network 
graph plays an important role in the definition of Thiemann's regularization of the Hamiltonian 
constraint [T3] . 

Both the Rovelli-Smolin and the Ashtekar-Lewandowski proposals, defined here on the node 
Hilbert space, admit classical versions: we dequantize the operators E\ to obtain vectors A\ S M 3 . 
This results in two distinct functions on phase space 

Ves(Ai) and V AL (A)- (4) 

Recently, a third proposal for the volume operator at a node has emerged [TTj- Motivated by 
the geometry of the Minkowski theorem, Bianchi, Dona and Speziale suggest the promotion of the 
classical volume of the polyhedron associated to {A{\ to an operator 

V po i y {Ai) -> V P oi y (Ei). (5) 

The number N of links at the node determines the number of faces of the polyhedron. One advan- 
tage of this proposal is that it is closer in structure to the spin foam formulation of the dynamics 
of loop gravity [181 CES] ■ 



In the case of a node with four links, N = 4, all three of these proposals for the volume operator 
coincide and match the operator introduced by Barbieri [19] for the volume of a quantum tetrahe- 
dron. The heart of this paper is a study of the semiclassics of this operator. Remarkably, the space 
of convex polyhedra with fixed face areas has a natural phase space structure. This was first shown 
by Kapovich and Millson [20J in a different context (see also |17[ [2T| [22] ) . Using this kinematics we 
study the classical volume dynamics and perform a Bohr-Sommerfeld quantization of its spectrum. 
The semiclassical analysis provides a useful synthesis and simplification of previous results: we find 
excellent agreement with the studies of the volume operator spectrum El [3 El 13 Ell E] and 
provide a new derivation of many of its features. Moreover, our analysis provides new geometrical 
insights on analytic properties of the volume spectrum, such as its degeneracy and the values of 
the maximum and minimum eigenvalues; for a detailed list of the new results see the conclusions 
in section [3 

We briefly recall the central ideas of the semiclassical methods used in this paper. The later 
development of Bohr's correspondence principle by Sommerfeld and Ehrenfest led to an elegant 
approximate quantization, now called Bohr-Sommerfeld quantization. The technique begins by 
finding the dynamical orbits of the system. In our case, the role of the Hamiltonian is played 
by the volume and we will analyze the volume orbits in phase spaceQ The Bohr-Sommerfeld 

1 More precisely, we use the oriented volume square Q to simplify the analysis. 
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quantization condition is then expressed in terms of the action integral I associated to each of 
these orbits: 

1(E) = jpdx = 2irk(n+ ~), (6) 

where the integral is over the orbit of energy E and the quantization level is denoted by n . Solving 
for E n gives a semiclassical approximation for the volume spectrum. 

It is interesting to notice that we are deriving the spectrum of the spatial volume in quantum 
gravity using methods that are older than quantum mechanics itself. In the 'old quantum theory', 
the quantization condition ^ was based on a physical assumption about robustness of discrete 
spectra under slow perturbations. Already several systems were known to have discrete equispaced 
spectra, for instance the harmonic oscillators in Planck's work on the black body spectrum, or in 
Einstein and Debye's work on specific heats. If the parameters of these systems are slowly modified, 
the energy spectrum changes but its discreteness is not destroyed. Lorentz and Einstein asked about 
a classical version of this robustness of spectra: i.e. whether there was a classical quantity that 
could be preserved under a slow — adiabatic — change of the external parameters of such systems? 
This question led Ehrenfest to identify adiabatic invariants as quantities that, on allowed orbits, 
are realized in multiples of 2irH. The Jacobi action integral 1(E) can be shown to be an adiabatic 
invariant for the system, and provides a natural candidate to extend Planck's quantization of the 
energy of the harmonic oscillator to other periodic systems, for example the Hydrogen atom. These 
methods continue to offer new insights into quantum physics by exposing the classical geometries 
underlying quantum phenomenon in the semiclassical limit. 

The structure of the paper is as follows: We begin by summarizing the setup for the node Hilbert 
spaces %n for general valency N. Next we describe how this space limits to a classical phase space 
and describe its Poisson structure. In section[2]we briefly describe what is known about the volume 
spectrum in loop gravity, focusing on the case N = 4 of the quantum tetrahedron. In sections [3] 
and [4] we describe the phase space of a classical tetrahedron, study the volume dynamics, and take 
up its Bohr-Sommerfeld quantization. The spectrum and its degeneracies are reported in sections 
[4] and [5} Elliptic functions play a prominent role in this quantization and section [6] develops tools 
for analyzing the largest and smallest volume eigenvalues using the properties of elliptic functions. 
Section [7] summarizes our results. 



2 Volume operator in loop gravity: the 4-valent node case 

The graph Hilbert space 7ir of loop gravity can be built out of gauge-invariant Hilbert spaces 
associated to nodes of the graph. For the entirety of this paper we focus on a single node n and 
its Hilbert space %n- Let us assume that the node is iV-valent and call the spins labeling its links 
jl (1 = 1,... ,N). To each representation ji we associate a vector space Hj t that carries the action 
of the (577(2) generators J\. The standard basis \jimi) is labelled by eigenvalues of the Casimir 
Jf = Ji ■ Ji and the z-component of J/ , J\ z . The space T~Ln is defined as the subspace of the tensor 
product 7ij 1 ® • • • (g) %j N that is invariant under global SU (2) transformations (the diagonal action) 

n N = Inv (% h ® • • • ® H jN ) . (7) 

We call %m the space of intertwiners and when emphasizing the parametric dependence on j\ , . . . , jn 
we write JLnUi, ■ ■ ■ , jjv)- The diagonal action is generated by the operator J, 



N 

j=Y^t 

i=i 



(8) 
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States of "Hat are called intertwiners and can be expanded in the \jimi) basis described above, for 
\i) G 7i N , 

\i) = Y^i mi - mN \hmi) ■■■\j N m N ). (9) 

m's 

The components i m ^'" m N transform as a tensor under SU(2) transformations in such a way that 
the condition 

J\i) = (10) 

is satisfied. These are precisely the invariant tensors captured graphically by spin networks. In 
loop gravity a special role is played by the operators 

Ei = 8ttj£ 2 p Ji, (11) 

as illustrated below0 

Let us now specialize to the case N = 4, a node with four links. The Hilbert space Ti^ is the 
intertwiner space of four representations of SU(2), 

n 4 = Inv (n h <g> U h ® H j3 ® H jA ) . (12) 

We introduce a basis into this Hilbert space using the recoupling channel 7ij 1 <g> Hj 2 and call these 
basis states \k). The basis vectors are defined as 

\ k ) = ^2 i t im2m3mA \h,mi)\j 2 ,m 2 )\h,rn ? ,)\j A ,m A ), (13) 
mi ■■■rriA 

where the tensor j™i m 2m.3"i4 j g defined i n terms of Wigner 3j-symbols by 



k 



■m 1 m 2 m3m4 
l k 



y (-if- m ( h 32 k )( k 33 k \ (14) 

^— ' V mi 7B2 m I y —m 7713 I 

m=—k 

The index k ranges from k m \ n to k max in integer steps with, 

kmin = max(|ji - j 2 \, \j 3 - j 4 |) and k ma _ x = min(ji + j 2 , h + k)- (15) 

The dimension d of the Hilbert space is finite and given byQ 

d = A; max — k m [ n + 1. (16) 

The states \k) form an orthonormal basis of eigenstates of the operator J% ■ J 2 . Following Barbieri 
[191 IHl E2] , these states can be understood as describing quantum tetrahedra. The operator E r ■ E s 
measures the dihedral angle between the faces r and s of the quantum tetrahedron |23j . 

The operator \J E T ■ E r measures the area of the rth face of the quantum tetrahedron and states 
in H4 are area eigenstates with eigenvalues 87rj£ 2 p ^ j r (j r + 1), 

E r • E r \i) = ^ZpsfjAjr + 1) (17) 



2 Here £p is the Planck length and 7 is the Barbero-Immirzi parameter. They should both be understood as 



coupling constants of the theory. Throughout the remainder of the paper we will take ip = 7 = H = 1 
3 This dimension can also be expressed in a symmetrical manner that treats all four j r on an 
d = min(2ji,2j 2 ,2j 3 ,2j4,ji + 31 + j 3 -jiji + h ~ h + H, h ~ h + h + .74, -ji +32 + jz + 3 a) + 1- 
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The volume operator introduced by Barbieri is 

V2 



V = ^-yJ\E 1 -(E 2 xE 3 )\, (18) 
and because of the closure relation 

(E X + E 2 + E 3 + E 4 )\i) = (19) 

this operator coincides with the Rovelli-Smolin operator for a = 2y/2/3. It also follows from ( |19[ ) 
that the Ashtekar-Lewandowski operator on H.^ is simply given by Val = vW\Vrs, where a is a 
number depending on the Grot-Rovelli class, see |24[|25j . of the link tangents at the node that can 
attain the values a = 0, ±1, ±2, ±3, ±4^] Therefore the Ashtekar-Lewandowski volume operator 
coincides numerically with the operators ^ and §5§ when the tangents to the links fall into the 
class corresponding to a = ±4 and the constant a is choosen to be 2\/2/3 as before. Otherwise, it 
is proportional to it. The volume operator introduced by Barbieri can be understood as a special 
case of the volume of a quantum polyhedron discussed in |17j . 

In order to compute the spectrum of the volume operator, it is useful to introduce the operator 
Q defined as 

Q = \e x ■ (E 2 x E 3 ). (21) 

It represents the square of the oriented volume. The matrix elements of this operator are easily 
computed and we report them momentarily. The eigenstates \q) of the operator Q, 

Q\q) = q\q), (22) 

are also eigenstates of the volume. The eigenvalues of the volume are simply given by the square- 
root of the modulus of q, 

V\q) = y/\q\\q). (23) 

The matrix elements of the operator Q in the basis | k) were originally computed independently 
by Chakrabarti and then by Levy-Leblond and Levy-Nahas |26[ [27] . They are 

here we introduce the shorthand A\ = ji + 1/2, indeed in the semiclassical limit the operator Ji 
with Casimir j\ is associated to an angular momentum vector A\ with magnitude A\ = ji + 1/2. In 
fact, as explained in [28J, semiclassical quantization of these angular momenta, which in this case 
is exact, gives a discrete area spectrum coincident with that of loop gravity. The function A(a, 6, c) 
returns the area of a triangle with sides of length (a, b, c) and is conveniently expressed in terms of 
Heron's formula 

A(a, b, c) = - yf{a + b + c)(a + b - c)(a - b + c)(-a + b + c). (25) 
4 The number a is defined as 

/ 1 2 3 4\ /12 3\ / 1 2 4\ . / 1 3 4s /23 4n , nn \ 

a(v ,v ,v ,v ) — e{v ,v ,v ) — s(v ,v ,v ) + e(v ,v ,v ) — e(v ,v ,v ). (20) 
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Computing the spectrum of Q amounts to computing the eigenvalues of a d x d matrix, where d is 



the dimension of the Hilbert space given in (16). This can be done numerically and several of our 



figures compare the eigenevalues calculated in this manner to the results of the Bohr-Sommerfeld 
quantization, see Sections [4] and [6j 

There are a number of properties of the spectrum of Q (and therefore of V) that can be 
determined analytically. We list some of them below and refer to Levy-Leblond and Levy-Nahas 
[27] for a detailed analysis: 

• The spectrum of Q is non-degenerate: it contains d distinct real eigenvalues. This is a 
consequence of the fact that the matrix elements of Q in the basis \k) determine a d x d 
Hermitian matrix of the form 



(26) 



/ iai • • • \ 
— ia\ iai ' ' ■ 
-ia 2 '•■ 

v : '•• '•• '•• / 

with real coefficients a^. For a derivation of this result see Appendix A. 

The non-vanishing eigenvalues of Q come in pairs ±q. A vanishing eigenvalue is present only 
when the dimension d of the intertwiner space is odd. These two properties are a direct 
consequence of the structure of the matrix (26) discussed above. This matrix is i times 



an antisymmetric matrix and the eigenvalues of antisymmetric matrices have this property. 
Physically, this is a consequence of the action of parity. 

From the above we see that the volume spectrum of a fixed node space Hn, that is, of an 
intertwiner space with given ji, ■ ■ ■ , j'4, is twice degenerate, 



V\±q) 



±q). 



(27) 



There is a second, exact degeneracy in the volume spectrum that has not been previously 
noted. Let s = l/2(ji + ji + j% + j'4) and j[ = s — ji (/ = 1, . . . , 4). The volume spectrum of 
the node spaces H^ji, . . . , j'4) and % 4 (ji, • • • , j 4 ) are identical. This is a manifestation of the 
Regge symmetry, which is briefly discussed in Section [5} 



• For given spins Ji, • • • , J4, Brunnemann and Thiemann have estimated the maximum volume 
eigenvalue using Gershgorin's circle theorem [9] and they find that it scales as v n 
where j max is the largest of the four spins ji . 



.3/2 

~ Jmaxi 



• They have also estimated the minimum non- vanishing eigenvalue (volume gap) and find that 

it Scales as V min ~ jmax [ll E] . 

Our semiclassical analysis reproduces all of these results and provides several new insights into 
their structure. For a parallel list of the semiclassical results see the conclusions. This completes 
our review and we turn now to the classical geometry of tetrahedra. 



3 Tetrahedral volume on shape space 

The finite dimensional intertwiner space Hn can be understood as the quantization of a classical 
phase space [TTJ [2TJ [22] . We set up this correspondence in three steps: first the relationship between 
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the algebra of vectors and the geometry of convex polyhedra (Minkowski's theorem) is explained, 
next we describe how to endow these vectors with Poisson and symplectic structures by interpreting 
them as angular momenta and, finally, we briefly describe how quantization of these objects leads 
IoTIn- The section concludes with the specialization of these results to the case of tetrahedra and 
a discussion of the classical volume. 

Minkowski's theorem, [29], states the following: given N vectors G R 3 (1 = 1,..., N) whose 
sum is zero 

A x + • ■ • + A N = (28) 

there exists, up to rotations and translations, a unique convex polyhedron with N faces associated 
to these vectors. Furthermore, the vectors A\ can be interpreted as the outward pointing normals 
to the polyhedron and their magnitudes A\ = \Ai\ as the face areas. Minkowski's proof is not 
constructive and so this is strictly an existence and uniqueness theorem. We call the process of 
building one of these polyhedra (given the area vectors) a Minkowski reconstruction. In the case of 
the tetrahedron the reconstruction is trivial but for larger N this is a difficult problem |17j . These 
convex polyhedra are our semiclassical interpretation of the loop gravity grains of space. 

Following Kapovitch and Millson [20j , we can associate a classical phase space to these polyhe- 
dra. We interpret the partial sums 



Pk 



k+1 

1=1 



(k = l,...,N-3) (29) 



as generators of rotations about the p k = A\ + • • • + A k+ \ axis. This follows naturally from 
interpreting each of the A[ vectors as a classical angular momentum. There is a well-known Poisson 
structure on a single angular momentum (the Lie-Poisson bracket, see [30]) and this extends to the 
bracket 

on N distinct angular momenta, here / and g are arbitrary functions of the A%. With this Poisson 
bracket the p k do, in fact, generate rotations about the axis p k = A± + - ■ - + A k +i. This geometrical 
interpretation suggests a natural conjugate coordinate, namely the angle of the rotation. Let q k be 
the angle between the vectors 

Vk=Pk* A k+1 and w k = Pk X A k+2 . (31) 
It is an easy check to show that 

{Pk,Ql} = hl- (32) 

The pairs (pk,qk) are canonical coordinates for a classical phase space of dimension 2(N — 3); 
this dimensionality is further explained below. We call this the space of shapes and denote it by 
V(Ai, . . . ,An) or more briefly Vjy. 

The explicit inclusion of the A\ in V(A\, . . . ,Ajy) highlights a parametric dependence of the 
space of shapes on the magnitudes A\. This parametric dependence is more natural when viewed 
from another perspective on V?q, i.e. viewing it as a symplectic reduction of the product of N two 
spheres, (S 2 ) N . Upon fixing the magnitude A\, the angular momentum vector A\ is restricted to a 
sphere of radius A[. This sphere is a symplectic leaf of the Poisson manifold described above and 
the collection of all the spheres (S 2 ) N can be endowed with the product symplectic structure. If 
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we now symplectically reduce (S ) by the zero level set of the momentum map 

N 



A = ^A l (33) 



i=i 

we once again obtain Vn- This also explains the dimension of the reduced space; because zero is a 
fixed point of the group action we lose twice as many dimensions as there are components of the 
momentum map, dim "Pat = 2iV — 6 = 2(N — 3). 

Finally, that the quantization of Vn is the Hilbert space Hn of an iV-valent node n can be 
seen as follows: Lift the N angular momentum vectors Ai via the Schwinger map to a phase space 
of 2iV harmonic oscillators, C 2N . These oscillators can be quantized in the standard fashion and 
once again the Schwinger map can be used, on the quantum side, to reduce to N quantum angular 
momenta. The constraint A\ + • • • + Afq = becomes 

Ji + • • • + J N = 0, (34) 

which is precisely the gauge invariance condition (81). 

We turn now to the classical analog of formula ( 181) , the volume of a tetrahedron as a function on 
the shape phase space V(A\, . . . , A4) = V4. This will be the starting point of our Bohr-Sommerfeld 
analysis in the next section. 

As discussed above the Minkowski theorem guarantees the existence and uniqueness of a tetra- 
hedron associated to any four vectors A\, (I = 1, . . . , 4) that satisfy A\ H V A& = 0. Without loss 

of generality we will take A\ < A2 < A3 < A4. In terms of these magnitudes, a condition for the 
existence of a tetrahedron is that A± + A2 + A3 > A4, equality yielding a flat (zero volume) tetra- 
hedron. This is clearly necessary, as there would be no way to satisfy closure if Ai + A2 + A3 < A4 
held. It is not difficult to argue that this is also sufficient for there to exist at least one tetrahedron 
with these face areas (in fact, there are infinitely many). The space of tetrahedra with four fixed 
face areas V(Ai, A2, A3, A4) = V4 is, as we will now argue, a sphere. 

Following the general construction outlined above, the canonical coordinates on V4 are p\ = 
\A\ + Ai\ and q±, the angle between v± = A\ x A2 and wi = (A± + A2) x ^3. For the remainder 
of the paper we adopt the simplified notation A = A\ + A2, A = \Ai + A2I = pi, and (j) = q\. 
Recalling that the Ai, (/ = 1, . . . , 4) vectors are to be thought of as generators of SU(2) actions we 
observe that A generates rotations of Ai and A2 about the ^4-axis. This action rotates vi about 
the, perpendicular, A-axis while leaving w\ fixed and thus increments the angle (p. This is the 
geometrical content of the Poisson bracket relation {A, (j)} = 1. Because A is fixed by this rotation 
the closure condition, 

A 1 + A 2 + 13 + 14 = 1+13 + 14 = (35) 

is also unaffected by such a rotation. Consequently, we have a whole circle of distinct tetrahedra 
for each value of A with A m i n < A < A max , where A m i n = max {A2 — Ai, A4 — A3} and A ma _ x = 
mm{A2 + Ai, A4 + A3}. The collection of these circles over the interval of allowed A values is the 
slicing of a sphere into lines of latitude over the range of its z diameter. This is made precise by 
the behavior at the ends of the range of A; either the vectors Ai and A2 or the vectors A3 and 
A4 become colinear, and consequently the four j4-vectors are coplanar. Such configurations are 
all equivalent up to overall rotations in M 3 and hence under the Kapovich-Millson reduction they 
correspond to a single point in the reduced space. As usual, the (p coordinate becomes ill defined 
at these poles, here this is because either v\ or wi vanishes. 

The shape space sphere is depicted in Figure [lj The embedding space is M 3 and can be thought 
of as a copy of the SU(2) Lie algebra associated to A, we will call this angular momentum space 
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Figure 1: The space of shapes for a tetrahedron, V(%, §, §, |). The darkened contours are quantized 
level sets of the classical volume squared, Q. One quantized level set is hidden from view. The 
inset depicts a tetrahedron corresponding to the dot on the q 3 level set. The point corresponding 
to the tetrahedron with the largest possible volume, given the face area constraints, is also marked 
with a dot. The shape space coordinates are indicated at the top of the sphere. 



A. Coordinates on A and the shape space sphere are defined as follows: place the origin of A 
at the center of the range of A and choose cartesian coordinates A x , A y and A z with A z = 
A — (A max + Amin)/2. The radius of the shape space sphere is R = (A max — A m \ n )/2 and the 
relation A z = RcosQ together with the angle eft from above define a spherical coordinate system 
on the sphere. The phase space geometry discussed here is a close analog to that of the symmetry 
reduced phase space of the 6j-symbol, see jH] for extensive discussion of the 6j-symbol case. 
The classical analog of the volume, (Il8|), is 



V = ^^/\A 1 .(A 2 xA 3 )\, (36) 

and just as with the quantum theory it will be more straightforward to work with the squared 
classical volume, 

Q = ~li ■ (A 2 x A 3 ). (37) 

A quick vector algebra calculation will convince you that these are indeed the volume and squared 
volume of a classical tetrahedron. Because Q is a rotational invariant quantity it projects onto the 
shape phase space and can be thought of as a function of the A and (f> coordinates, Q(A, (j))^ This 
expression is easily derived by computing v± x wy, from the definitions of v\ and w\ one finds 

v x xwx = ^QA. (38) 
Note that the magnitude \vi\ = \A\ x A 2 \ is equal to twice the area A of a triangle with side lengths 



A\,A 2 and A and, using the closure relation (35), similarly = |^3 x A±\ is twice the area A 



5 Because A and A z differ by a constant shift we will freely switch between them. 
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of a triangle with side lengths A^,A^ and A. The definition of 4> as the angle between v\ and w\ 
allows us to conclude that the volume squared is 

Q=^sin0. (39) 



Calculating the areas A and A using Heron's formula (25), 



A = \ ^/[(A 1 +A 2 ) 2 -A 2 ][A 2 -(A l -A 2 ) 2 ], (40) 
A = \yJ[{M + A,) 2 - A 2 ][A 2 - (A 3 - M?\, (41) 
shows that for fixed Ai,...,Ai, Q is indeed only a function of the coordinates A and 4>. The 



expression (39) will be the central tool of our Bohr-Sommerfeld quantization. In anticipation of 



the results of the next section, some quantized level sets of Q are depicted in Figure [Tj 

4 Bohr-Sommerfeld quantization of tetrahedra 

As discussed in the previous section, the phase space of a tetrahedron is two-dimensional with 
canonical coordinates A, cp, 

{A,4>} = 1. (42) 



The oriented volume square Q(A, <f>), Eq. (|39j), is a function on phase space and generates an 
Hamiltonian dynamics in a parameter time A. The Hamilton equations are given by the familiar 
formulae 

S = M • fx ={*«■ < 43 > 

Several of the orbits generated by Q have been displayed in Figure [TJ Along these orbits Q is 
preserved 

Q(A, <j)) = q. (44) 

As discussed in the introduction, the Bohr-Sommerfeld quantization condition is expressed in terms 
of the Jacobi action integral I associated to each of these orbits: 



I(q) = j> Ad4> = 2ir(n+-), (45) 

here we have denoted the level value of Q by q, the quantization level by n, and we take units in 
which h = 1. In the main body of the section we calculate I and impose this quantization condition. 
This is the usual Bohr-Sommerfeld procedure with Q playing the role of the Hamiltonian H and q 
the role of the energy E. 

We commence by solving the volume dynamics of A. The classical evolution of the tetrahedron 
with Q taken to be the Hamiltonian turns out to be easiest to calculate in terms of A 2 , 

^ = {A 2 , Q} = 2A{A, Q} = 2A^, (46) 

where A is defined to be the variable conjugate to Q and <\> is, as in the previous section, the angle 



conjugate to A. The right hand side can be evaluated by differentiating (39) with respect to 

d(A 2 ) 16 



d\ 9 



AAcos</>. (47) 
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Equation (39) can be used again to eliminate the cosine function, 

d{A 2 ) 1 



dX 



(4A) 2 (4A) 2 - (2A) 2 (9Q) 2 . (48) 



The argument of the square root plays an important role in what follows and so we introduce a 
shorthand for it, 

P{A 2 , Q 2 ) = (4A) 2 (4A) 2 - (2A) 2 (9Q) 2 . (49) 

This is a quartic polynomial in A 2 and while the general expressions for its roots are complicated 
they simplify when Q = 0, and so we further define, 

P(A 2 , Q 2 ) = P(A 2 , 0) - {2A) 2 {9Q) 2 = P (A 2 ) - (2A) 2 (9Q) 2 . (50) 

Equations ( |40[ ) and (41) yield a factored expression for Pq(A 2 ): 

P (A 2 ) = [A 2 - (A, - A 2 ) 2 ][A 2 - (A 3 - A 4 f][{A x + A 2 ) 2 - A 2 ][(A 3 + A 4 ) 2 - A 2 }. (51) 



Taking the shorthand x = A we can separate variables in (48) and integrate to find, 

X(x) = 9 dX =, (52) 

Jr 2 V (x - n)(x - r 2 )(r 3 - x)(r 4 - x) 

we assume that the four distinct real roots (ri , r 2 , r 3 , r 4 ) of the quartic P(x = A 2 , Q 2 ) are ordered 
as t\ < r 2 < T3 < r4. This is an elliptic integral; to bring it to the standard Jacobi form we use a 
Mobius transformation that brings the quartic to a conventional one with roots ±1, ±^=|^| This is 
possible as long as the cross-ratio of the r's is the same as the cross-ratio of the conventional roots, 
we use this to set the elliptic parameter m. Explicitly the substitution is 

z 2 = {U ~ ^ " X \ and m= ^-r2){r,-r l ) 
{rs - r 2 ){r 4 - x) (r 4 - r 2 )(r 3 - ri) 

After evaluation of the integral and some algebra this leads to the solution, 

X(A) = A 2 (X) = ^-^-r^-r^-^m)^ 
(r 4 - r 2 ) - (r 3 - r 2 )sn 2 (^, m) 

with 

V( r 4 - r2){r 3 - n) 

This is a complete solution of the dynamics. After the specification of a value for the volume, 
the quartic P(x,Q 2 ) can be solved and the volume evolution of the intermediate coupling A is 
given by (54). The evolution is periodic and the period can be expressed in terms of the complete 
elliptic integral of the first kind K(m) by T = 9g x 2K = 18gK. The fundamental period of 
the elliptic functions, To = 4K, is halved because they appear squared. For definiteness, in what 
follows we will assume that the elliptic parameter m is less than one. If this is not the case apply 
the transformation sn(u,m) = ^^sn(^/rnu, ^=) and you will find that the effect on (54) and (55) 
is to switch the roles of n and r 2 throughout. 

6 To avoid notational conflicts with the intertwiner eigenstates \k) all elliptic functions are written in terms of the 
elliptic parameter m = k 2 instead of the elliptic modulus k. 
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Before proceeding to the calculation of the action I of a curve A(X), we pause to describe some 
of the properties of the quartic P(x, Q 2 ) that will be useful in this calculation. For the value Q 2 = 



the quartic simplifies and is given by Pq(x). In particular Pq{x) can be explicitly factored, see (51 ), 
and we will call its four positive real roots f\ < f 2 < f$ < f±. For non-zero real values of the 
volume Q, the roots of the quartic equation P(x, Q 2 ) = are given by the intersections of the line 
y = 32AQ 2 x with the quartic y = P(x,0) = Po(x). This leads to a few general remarks about the 
roots {n, r 2 , ?"3, of P(x, Q 2 ). The largest and smallest roots are always real and satisfy < r\ < 
fi = min{(Ai -A 2 ) 2 , (A3-A4) 2 } and r 4 > f 4 = max{(Ai + A 2 ) 2 , (A 3 + A 4 ) 2 }. Meanwhile for small 
enough Q the middle two roots are also real and satisfy, r 2 > f 2 = max{(Ai — A2) 2 , (A3 — A4) 2 } 
and r% < f% = m\n{(Ai + A2) 2 , (A% + A4) 2 }. As Q grows the middle two roots coalesce and then go 
off into the complex plane. These observations are summarized graphically in the first two panels 
of Figure [5} 

The roots of a polynomial coalesce when the polynomial and its first derivative simultaneously 
vanish, P(x) = P'(x) = 0. With the notation introduced above P(x) = Pq(x)—?>2AxQ 2 (we suppress 
the Q dependence), two roots will coales ce w hen P(x) = and dP/dx = = Pq(x) — 324Q 2 or 



Qcoal = 1/18a/Pq(x) both hold. In section 6.1 we will show that the latter is precisely the condition 



needed to achieve the maximum volume of a tetrahedron with four fixed face areas, 



Qmax = I/I^P^X). (56) 

This means that the quartic roots coalesce precisely when the maximum real volume of the tetra- 
hedron is achieved and so for real volumes we need only consider real positive roots. Below we find 
that the action is given by complete elliptic integrals and it will be useful to be able to assume that 
the roots that arise in these formulas are real and positive. 

The action integral can now be calculated. The action for an orbit 7, e.g. one of the curves of 
Figure [TJ can be re-expressed in terms of the conjugate variable A to the volume, 



I = <p A{(j))d(p = <p A(X)-^dX. (57) 



Once again turning to (39) and solving for <j) gives, (f> = arcsin (9^4Q/(8AA)) and differentiating 
with respect to A yields, 



#_ 1 / 9Q dA 9AQ d(AA)\ 

d\ ~ y/i - (9AQ/8AAy \8AAdX ~ 8(AA) 2 dX ) 



(Q _ Q d{AA) \ 
\A AA dA J ' 



(58) 



Returning to the expressions for A and A ((40) and (41)) one can calculate d(AA)/dA and obtain, 



A 2 A 
A^-dX = <b QdX - (b Q [ — — — rTpr + 



dX J ^ J^y(A 2 -(A 1+ A 2 ) 2 ) ' (A 2 - (A, -A 2 ) 2 ) 

A 2 A 2 \ 

+ {A 2 _ (A3 + Ai?) + {A 2 _ {As _ A4 )2) ) A - ( 59 ) 

Because Q is a constant along the orbit, the first integral simply gives the period of the elliptic 
function, 

(j> QdX = 18gQK(m), (60) 
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where again the elliptic parameter m is given by (53). The remaining integrals are all of the 
same type, we introduce a parameter fj which takes the values of the roots at zero volume, (A\ — 
A 2 ) 2 , (A 3 - A 4 ) 2 , (A 1 + A 2 ) 2 , (A 3 + At) 2 , respectively for i = 1, ... ,4, and find, 



r3(r 4 - r 2 ) - r 4 (r 3 - r 2 )sn 2 (^,m) 
(r3 - f;)(r 4 ~ r 2 ) - (r 4 - n)(r 3 - r 2 )sn 2 (^,m) 



(61) 



These integrals are more simply expressed in terms of u = X/9g and the integration over a complete 
period is over the interval u 6 [0, 2K\. We find, 



n A 2 18 5 r 3 
Q (A^T) dA = - 



1 



(?"3 - Jo 1 - afsn 2 (u, m) 



du 



Q 



lSgr^af 



A 



sn iu,m) 



where 



a,- 



(r 4 - n) Jo 1 - a 2 sn 2 (u, m) 
(r 4 - fj)(r 3 - r 2 ) 



du, 



(r3 - rj)(r 4 - r 2 ) 
Collecting all four of these integrals we have, 

i 



(62) 



(63) 



1 



(?"3 - r i) Jo 1 - a 2 sn 2 (u, m) 
K sn 2 (u, m) 



du 



r 4 a 2 



(64) 



(r 4 - rA y 1 - a 2 sn 2 (u, m) 
which can be evaluated in terms of complete elliptic integrals yielding, 



du 



lSgQ 



'•4 



n(X4 - r 3 ) 
■[ (r4 - fj)(f3 - n) 



n(a 2 , m) 



(65) 



where n(a 2 , m) is the complete elliptic integral of the third kind. To highlight the structure of this 
result we condense the dependencies on the roots into two coefficients, 



a = 18s 



1 (»"4 



and 6j = 



18grj(r 4 - r 3 ) 
(r 4 - f;)(r 3 - fi) 



(i = l,...,4). 



(66) 



This allows us to write / in the more compact form, 



I = ^aK(m) - ^n(a 2 , m) \ Q. 



(67) 



This is our main result. Figure [2] displays a plot of this function for the same parameters used in 
Figure [T| 

In spite of the complex manipulations used to obtain the expression (67) its interpretation is 
simple: Applying Stoke's theorem to the action / = <f Ad(p of an orbit 7, we can interpret this 
integral as the symplectic area contained within the orbit. The symplectic form on shape space is 
determined by the Poisson bracket relation {A, (f\ = 1 and is oj = dA A d(j) so that 



A(cj))d<p = uj 



dA A 1 



(68) 
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Figure 2: A plot of the action integral I(Q) for the same parameters as Figure [I] The quantized 
levels q n shown satisfy the Bohr-Sommerfeld quantization condition. The corresponding orbits are 
shown in Figure [TJ 



We can just as well work with dA z A dip = Rd cos 9 A dip because A and A z differ by a constant and 
so this symplectic area only differs from the solid angle on the sphere by a normalization factor R. 

Although this action has a closed analytic form, finding the Bohr-Sommerfeld spectrum requires 
a numerical inversion. Recall that the strategy is to find the volumes for which the corresponding 



orbits capture (n + l/2)2-7r worth of area on the sphere. The analytic expression (64) has a compli- 
cated dependence on the volume; it appears explicitly as an overall multiplicative factor but also 
implicitly through the elliptic function's dependence on the roots ri, . . . ,r&, all of which depend 
on Q. Instead of trying to invert the action analytically we have calculated its value for several 
hundred points in the range of classically allowed volumes, interpolated between these values and 
numerically found the volumes for which the Bohr-Sommerfeld condition is satisfied, this process 
is illustrated schematically in Figure [2} The results of this analysis are presented for two examples 
in Figure [3] along with the numerical diagonalization of the volume matrix elements discussed in 
secti on [2] Recall that V = \/\Q\ and that the eigenvalues satisfy the same relationship v = ^/\q\ 
(see J23p. 

The Bohr-Sommerfeld approximations developed above effectively reproduce the volume spec- 
trum of loop gravity, see Figure [3j In a manner that is characteristic of this type of semiclassical 
approximation, the agreement rapidly improves with increasing quantum numbers, in our case 
increasing j[. This is illustrated in Appendix [Cj which tabulates the Bohr-Sommerfeld and loop 
gravity volume eigenvalues for a variety of spins. Section [6] contains additional comparisons of the 
numerical and Bohr-Sommerfeld results. The reduced quality of the approximation for the case 
where all the spins j r are equal, exhibited in the right plot of Figure [3] at j = 4, is also discussed. 
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Figure 3: A comparison of the Bohr-Sommerfeld and loop gravity Volume spectra. On the left: 
configuration with spins {j, j, j, j + 1}. On the right: configuration with spins {4, 4, 4, j} and 
j varying in its allowed range. The Bohr-Sommerfeld values of the volume of a tetrahedron are 
represented as dots, the eigenvalues of the loop-gravity volume operator as circles. Recall that the 
spins and areas are related by A\ = ji + 1/2. 



5 Degeneracy of the volume spectrum and the Regge symmetries 

Besides accurately reproducing the frequency of many emission lines in atomic spectra, the old Bohr- 
Sommerfeld quantization is also able to predict the intensity of such lines, that is, the degeneracy 
of the energy levels. We now apply these ideas to the case of the volume spectrum. 

The volume dynamics generates closed orbits in phase space. In fact, for every allowed non-zero 
value of the volume there are two closed orbits satisfying V(A,cp) = v. The two orbits are sent 
one into the other by a parity transformation of the tetrahedron and while the oriented volume 
square Q(A, (f>) = q is odd under parity, the volume eigenvalue v = y/\q\ is invariant. Thus the two 
classical orbits with the same volume give rise to twice degenerate volume levels. The orbit with 
vanishing volume is invariant under parity and so, v = 0, when it is an allowed quantum level, is 
non-degenerate. This is consistent with the following observations about the classical phase space: 
The total symplectic area of our phase space is finite. The Bohr-Sommerfeld condition implies that 
this area is divided up into a finite, integer number of Planck cells, 



dA A . 



2-nd, 



(69) 



here the integer d is also the dimension of the Hilbert space, (16). When d is odd, the Bohr- 
Sommerfeld condition predicts that there are d = 2n + 1 levels: n doublets of non- vanishing volume 
and one singlet of vanishing volume. Figure [2] provides an example of this case. 

There is a second more subtle degeneracy: the tabulation of Appendix [C] exposes a degeneracy 
in the volume spectra of distinct intertwiner spaces. The clearest example of this in our tables, is 
for the spins (juh,j3,u) = (6,6,6,7) and j' 2 , j' 3 , = (^,^,^,^); the volume spectra of 
these two intertwiner spaces agree exactly. 
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This degeneracy can be understood in terms of a classical geometric version of the well known 
Regge symmetries of the 6j-symbol [32] • Let s = l/2{j\ +J2 + 33 + ji), if the spins of a 6j-symbol, 

{{ n H (70) 

[J3 J4 J23,J 

are transformed to j[ = s — ji, (I = 1, ... ,4) and j' 12 = ji2, j'23 = 323, then Regge showed that, 



{* 


32 


J12I 




J2 




{33 


k 


j'23/ 


u 


3\ 


323) 



(71) 



Roberts explains that this symmetry can be understood geometrically as a scissors congruence 
of the relevant tetrahedra [33]. Two polyhedra are scissors congruent if the first polyhedron can 
be sliced into finitely many polyhedral pieces and then reassembled into the second polyhedron. 
Evidently, scissors congruence preserves volumes. 

In the present context this symmetry completely explains the degeneracy of the volume spectra; 
besides the Minkowski tetrahedron having the vectors Ai as face normals, we prove this using an 
auxiliary tetrahedron having edge vectors A\. This auxiliary tetrahedron has volume 7^731 Q = 
^\ A\ ■ (A2 x A3). Now, performing a Regge transformation we obtain a second auxiliary tetrahedron 
that is scissor congruent to the first one. The two auxiliary tetrahedra have the same volume, 
Q! = Q, and this implies that the volumes of the corresponding Minkowski tetrahedra are also 
equal V = y/\Q'\ = \J\Q\ = V . Moreover, for every auxiliary tetrahedron in the phase space V4 
we have a scissor congruent tetrahedron in the Regge transformed phase space V 4 . As a result the 
action integrals also coincide, I(q) = I'(q), and the degeneracy of the volume spectrum is explained. 

The Regge degeneracy of the spectrum leads to a conjecture: We have shown above that two 
Minkowski tetrahedra related by a Regge symmetry have the same volume. Are they also scissor 
congruent? We conjecture that they are. The proof of this conjecture hinges on showing that 
these two tetrahedra have the same Dehn or Hadwiger invariants, again see [33]. We leave the 
investigation of this conjecture open for future work. 



6 Limiting cases: largest and smallest volumes 

6.1 Classical analysis of limiting cases 

Before proceeding to a Bohr-Sommerfeld analysis of the limiting values of the volume spectrum, 
we investigate the classical extrema of \Q\. At the classical level, the minimum of \Q\ is always 
zero: Open the angle <j) until the plane spanned by A\ and A2 coincides with the plane spanned by 
A3 and A4, then Q = 0. This, however, does not lead to the conclusion that the minimum of V is 
always zero. This is because the Minkowki theorem does not hold for planar configurations of the 
Ai, instead it is generically singular for these configurations. Certainly Q = implies that V = 0, 
however, the issue is the correspondence between a planar set of A\ and such a flat tetrahedron. 
Geometrically this is clear, a flat tetrahedron has faces that lie in a plane and the normals to these 
faces are all collinear. Thus it is only the subset of planar configurations of vectors Ai that are 
actually collinear that can have a Minkowki type correspondence with a flat tetrahedron. 

We can say more: the collinear vectors must satisfy closure and so, for some choice of signs we 
must have ±Ai±A2±A3±A4 = 0. Taking into account the ordering convention A\ < A2 < A3 < A4 
we can bring the number of cases down to just two 



A 2 - Ax = A A - A 3 or A 2 + A 1 = A 4 - A 3 . 



(72) 



6 LIMITING CASES: LARGEST AND SMALLEST VOLUMES 



17 



The latter condition leads to a trivial shape space consisting of a single point because A m i n = A max . 
We will call these two conditions the "flatness" conditions. They are also significant for the Bohr- 
Sommerfeld quantization. 

Note that, even for collinear configurations of the A\ the Minkowski theorem still doesn't hold. 
The trouble is uniqueness. An infinite number of flat configurations all share the same area vectors. 
In fact, the differential structure of the shape space breaks down when the flatness conditions 
are satisfied (for an analogous observation see |34|). The qualitative picture is that when flat 
configurations are present the phase space sphere develops a cusp and looks more like an inverted 
rain drop. 

Notice that everything that has been said up to this point is in regards to a tetrahedron 
that is exactly flat. This is significant because it highlights the singular nature of the Minkowski 
construction for planar configurations of the A\. One can construct a tetrahedron with arbitrarily 
small volume from a set of vectors that is arbitrarily close to planar, there is only trouble when 
exactly flat tetrahedra are desired. The precise treatment of flat configurations warrants further 
investigation. 

For the purposes of the present work we summarize the preceding observations: unless the 
flatness conditions are satisfied, Q = should not lead to the conclusion that there is a constructible 
tetrahedron with V = 0; if a flatness condition is satisfied then only the corresponding pole of the 
shape space sphere corresponds to flat configurations with V = and it generally corresponds to a 
whole class of such tetrahedrsj^j thus the great circle on which Q = (<ft £ {0) 71 "}) wm be regarded 
as mathematically useful but largely physically meaningless. 

Turning to the maxima of Q, take three faces to have fixed areas, say A\,A 2 and A3 but the 
full vectors A s , (s = 1, 2, 3) not to be given. Writing the triple product of Q as the determinant of 
a matrix M = (A±, A2, A3) whose columns are the vectors A%, A2 and ^3 and squaring yields, 

/ A\ a\-a\ A x -A 3 

Aj-Aj A 2 -A 3 \, (73) 

\A3-Ax A 3 -A 2 A\ 



Q 2 = — det M T det M = — det 
^ 81 81 



where M T denotes the transpose of M. Taking the unknown dot products A s ■ A t , (s < t = 1, 2, 3) 
as variables and extremizing one finds the minima already discussed, where A%, A2 and ^3 are 
collinear, and a single global maximum where A s ■ At = 0, (s < t = 1, 2, 3). This maximum must 
satisfy closure and so, 

A\ = A\ + A\ + A\ + 2Ai • A 2 + 2A X ■ A 3 + 2A 2 ■ A 3 = A\ + A\ + A\. (74) 

Then the maximum volume of a tetrahedron with three fixed face areas is the one with three right 
dihedral angles and the fourth face area given by the equation above. If the same technique is 
used to maximize the volume over the space where all four face areas are given then the closure 
condition must be implemented as a constraint. For fixed A\, A2 and ^.3 it is clear that the 
constrained maximum will not be larger than the one just found and will be equal to it when 



the fixed value of A± is that of (74). Rather than implementing the constraint with a Lagrange 



multiplier it is easier to extremize equation (39), which treats the A[,(l = 1, •••4) on an equal 
footing. Once again we use the shorthand Po(A 2 ) for the polynomial part of the area product A A, 

P (A 2 ) = [A 2 - (A! - A 2 ) 2 ][A 2 - (A 3 - A 4 ) 2 ][(A! + A 2 ) 2 - A 2 ][(A 3 + A,) 2 - A 2 }, (75) 



7 A simple example: Consider the flat tetrahedron with four equal areas Ai — |a 2 (/ = 1, ... ,4). One example 
would look like a square made out of two right triangles of side length a and with, say, upward pointing normals and 
two downward normals. However, this can also be achieved with any rhombus that has side length b and acute angle 
/3 satisfying b 2 sin /3 = a 2 . Thus there is a one parameter family of tetrahedra with Ai = |a 2 . 
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that is, A A = 1 /16 V /P (A 2 ). The expression for the squared volume ( |39| ) simplifies to Q 
l/(18y / x) \/Pq{x) sine/;. This expression is maximized when <j> = 7r/2 and when dQ/dA = or, 



P (x) = xP' {x). (76) 

This condition is another quartic equation, the roots of which are not worth explicitly displaying in 
general, however choosing the root which maximizes Q, say x, then we have Q max = 1/18-v/Pq(^)) 



which justifies our claim at (56). This expression will be useful below for finding the largest 
eigenvalues of V. These two cases exhaust the most natural constraints on the face areas. If you 
fix only two of the face areas the volume can grow without bound. 

So far we have considered only constraints on the face areas. Another natural constraint is to 
require that the total surface area of the tetrahedron be constant and look for the largest volume 
within this class. As one might expect, similar arguments to those presented above lead to the 
conclusion that the largest volume tetrahedron under this constraint is the one with all face areas 
equal. This concludes our general treatment of the extrema of Q. 

6.2 Largest eigenvalues 

The limiting behavior of the volume spectrum for large and small eigenvalues can now be explored 



with the assistance of the analytical formula (67). Let us first consider large eigenvalues: the 



volume function attains a maximum on the sphere and so our strategy will be to expand the action 
function around this maximum. We have, 



dS 

5(g) = S(q max ) + — 



(q - <?max) H • (77) 

Q = <?max 



From the theory of action-angle variables, the derivative of the action with respect to the Hamil- 
tonian (in this case the volume) is the period T of the system. For the largest eigenvalue, we are 
in the same situation as if we were finding a ground state, that is, instead of capturing 2ir worth 
of area on the sphere this state only captures an area tt and so, 

, s _ fffanax) ~ 5(g) _ 7T 

(gmax - q) - 7fT ( s — =7 s- l'°J 



Then the largest eigenvalue is given by, 



(gmax- ^J , )■ (79) 



At equation (76) we found that the maximum classical volume is attained when, 

P (x) = xP^{x) (80) 

or more explicitly when, 

1 _ 1 1 1 1 

x ~ x - (Ax - A 2 f + x - (A 3 - Atf + x - {A l + A 2 ) 2 + x - (A 3 + Atf ' (81) 

The roots of this quartic for generic A\ are complicated functions of the A\ , however in the case A\ = 
• • • = A4 = Aq this equation is easily solved and one finds, x = \A%. The maximum classical volume 
in this case is Qmax = (2 3 /3 7 / 2 )Ajj and the period is T(q mSLX ) = 18gK(0) = 18(v / 3/(2A 2 ))(vr/2) = 
3 5 / 2 7t/(2^4q). The maximum eigenvalue is given by putting these values into (79), 



- w^f^k- < 82 > 
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3 /2 

This reproduces the A Q scaling that has been found in previous works and refines it to the next 
order. This scaling is plotted as the uppermost line in Figure [7j Further corrections could be 
developed by retaining more terms in (77). As discussed in Section 6.1, this is the appropriate 



scaling to consider under the constraint of fixed total surface area. Consequently, this scaling may 
be interesting to investigate in more detail within the U(N) framework |35j . 

Occasionally the equal area tetrahedron has been supposed to be the one whose volume grows 
most rapidly as the areas are increased. However, we have seen above that this is not always the 
case. The tetrahedron with maximum volume depends on the space under consideration. For the 
space with three face areas fixed, corner tetrahedra with, for example, A\ = A<i = A% = Aq and 
A^ = \/3Aq are maximizing. Indeed, solving (79 ) for these corner tetrahedra we find a larger scaling 
coefficient, 



,V2 




(83) 



In Figure [4] the Bohr-Sommerfeld spectrum is compared to the numerically calculated exact spec- 
trum and the scaling derived here for the corner tetrahedra. 
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Figure 4: This figure compares the semiclassical scaling of the largest volume eigenvalue (dark line) 
to the Bohr-Sommerfeld (dots) and loop gravity (circles) spectra. The first three spins are given 
by ji = j'2 = j3 = j and and the largest spin j'4 is given by the closest integer or half integer to 
\/3j, that respects the Clebsch-Gordan conditions. 



6.3 Smallest eigenvalues: non-flat configurations 

The small eigenvalue case has subtleties associated with it. As is clear from the Taylor expansion 

dl 



I(q) = 1(0) 



0q 



■ + 



(84) 



q=0 



6 LIMITING CASES: LARGEST AND SMALLEST VOLUMES 



20 



the smallest eigenvalues are associated with the longest period at zero volume, q mm = (1(0) — 
I(q))/T(0). This would be the end of the story except for the fact that there are a number of shape 
spaces for which the period at zero volume can become infinitely long and the Taylor expansion 
above is invalidated. For this reason, we have to treat different shape spaces individually. We will 
begin by treating the case where the period is finite, as it is simpler. 

One subtlety of the small volume cases is immediate: the smallest eigenvalue depends on the 
dimension of the intertwiner space, d = dim ^4. As already noted, the volume is odd under parity 
and so there is only a zero eigenvalue when d is odd and this is the only volume state invariant 
under parity. We are interested in the first non-zero eigenvalue and thus when d is odd the spacing 
1(0) — I(q) is the spacing between two quantized orbits, h, or in our units 2ir. For the finite period 



case that we are considering the Taylor expansion (84) is valid and we have, T(0)g mm = 2ir. On 
the other hand if the phase space is even dimensional, neighboring quantized orbits evenly straddle 
the zero volume contour and we have 1(0) — I(q) = tt, so that T(0)q m i n = n. 

We found the period above, in general it is T = 18gK(m), where g and m are given in equations 
([53|) and (54). For q = the quartic P(x, Q = 0) factorizes and the roots are f\ = (A\ — A2) 2 , r% = 



(A3 — A^fyrz = (Ai + A2) 2 \r 4 = (A3 + A4) 2 , for example. Of course, depending on the choice of 
Ai, . . . , A4 other orderings are possible (recall that the r are defined such that r% < r<i < r% < r 4 ). 



If, as suggested after equation ( 54 ) , m is always chosen such that m < 1 these other orderings lead 



to the same result: g = 1 / (2>/ A1A2A3A4) and 

= (M + A2 + A3- A4XA1 +A2-A3 + A 4 ) 
(2A,)(2A 2 ) 

< (^1 -A2 + A3 + A 4 )(-A x +A2 + A3 + A 4 ) 
(2A 3 )(2A 4 ) 

The complete elliptic integral of the first kind has the power series expansion, 



K(m) = £ £ 



2 



(2n)\ n:! 
2 2n n! 2 



m n , (86) 



valid for m < 1. If m is small enough such that higher order terms can reasonably be neglected 
then, for d odd, 



and this can be improved as much as desired by including more terms from (86). Expressed in 
terms of the volume of the tetrahedron we have, 

, . , , . ,1 1 a I 2/3 if d is odd „ . . 

v mia « (A 1 A 2 A 3 A 4 ) 1 / 4 { ' . for m « 1. (88) 

I V2/3 it d is even, 

6.4 Smallest eigenvalues: flat configurations 

The exact volume eigenvalues derived by Brunnemann and Thiemann [9] are for special cases of 
the Ai,'" ,A 4 , such as A\ = A2 = 1 and A3 = A 4 = j + 1/2, where only two of the four vectors 
grow as you increase j. For cases like these we find the same qualitative scaling from the formula 
(88), v ~ J 1 / 2 . However, notice that this special choice of the A r leads to m = 1 and invalidates 



the expansion of the elliptic integral K(m), in particular K logarithmically diverges (hence also 
the period) and a different approach to estimating the eigenvalues is necessary. Brunnemann and 
Thiemann were led to consider these special cases by their numerics, they found that these were 
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the phase spaces which lead to the smallest overall values for the volume spectrum. Below we 
describe what is special about the geometry of these cases and develop an alternative technique for 
estimating the spectra of these spaces. 

The longest periods are achieved when m = 1 and the elliptic function theory limits to ele- 
mentary functions, for example, the standard roots of the Jacobi form of the elliptic functions, ±1 
and zkl/y/m, degenerate. We worked out the elliptic modulus in the zero volume limit above, see 



(85). Setting this expression equal to 1 we find four roots, which can be regarded as expressing 
any one of the A r in terms of the other three, these are: A\ — Ai = A^ — A3, A\ + A2 = A3 + A4, 
A\ + A2 = —A3 — A4 and A2 — A\ = A4 — A3, only the last of which is physical due to our order- 
ing convention on the A\. This final condition is the first of the flatness conditions introduced in 



Section 6.1 (The other flatness condition arises when you consider the m = limit.) Remarkably, 
the elliptic function theory limits to elementary functions precisely when the shape space contains 
flat configurations. This provides a geometrical interpretation to the special cases that Brunneman 
and Thiemann investigated^ [9] : they are special because of contributions from flat tetrahedra. 

This flatness condition also has a nice interpretation in terms of the polynomial Po(x); it 
indicates that two of the polynomials roots are coalescing. Put this together with our observation 
that the maximum volume is achieved when r<i = T3 and, consequently, when m = and we find 
that Pq{x) is quite useful for characterizing the shape space of a tetrahedron with fixed face areas. 
These findings are summarized in Figure [5j 

In the case where m = 1 we can no longer use the Taylor expansion and elliptic function period 
to find the eigenvalues, instead we have to find the small volume behavior of the action function and 
try to invert it. The action depends on the volume in two ways, an explicit dependence through the 



prefactor Q in (64), and an implicit dependence through the roots of the quartic P(x, Q). Thus the 
first step in finding the small volume behavior is to expand the roots as a power series in Q. This 
would be quite laborious if we had to go through the solutions to the quartic equation, happily a 
simple alternative exists because we know the roots f j at Q = 0. We simply plug fj + p into P(x, Q) 
and require that p is such that the equation is satisfied at lowest order in Q, this process can be 
iterated to find rj to the desired order in Q. The result of these calculations to fourth order in Q 
are summarized in Appendix B. 

Using the series expansions of the roots we can Taylor expand the action as a power series 
in Q. This is slightly delicate because the complete elliptic integrals diverge logarithmically at 
m = 1, however, taking care to expand the logarithms to the proper order we find simple results. 
First we consider the case in which the A\ are all equal with A\ = Aq, so that f 1 = f2 = and 
f 3 = f4 = f = (2Aq) 2 . The action simplifies and the expansion yields, 



I = 18gQ 



1 - 2- 



K{m) - 2 , rir t ( r3) - U(al m) ) (89) 
[7-4 — r)[r3 — r) 



[r4 — r) 

=s^ + ^ln('J^) 3 + 0(Q 2 (ln(3)) (90) 

='(°)+(6- i/2 >(j^Hjy. 0») 

where in the last equality we have recognized ^ffn = 2ttA as half of the symplectic area of the 
sphere and hence the action at zero volume squared, 1(0). 

At lowest order then, the relation between / and Q can be inverted using Lambert's W function. 
Because this inverse function is uncommon, we briefly review the properties used in this work. The 



8 The flatness condition explains most of their cases. A few cases are special and not due to the flatness condition 
but rather because the dimension of the space of intertwiners is two. 




Figure 5: This figure summarizes our characterization of the shape space using the quartics Po(x) 
defined by equations (75) and (50). The straight lines are given by y = 324Q 2 x and the 
P(x,0) 



and P(x 

quartic curves by y 



Pq(x). The upper left panel shows the generic case in which the 
fixed volume determines four distinct roots. The upper right panel shows the coalescence of the 
middle roots at maximum volume and m = 0. The lower two panels show the coalescence of other 
pairs of roots at zero volume and m = 1 (these are two distinct cases) . 



Lambert W is defined as the function which inverts the relationship, 

x = We w , (92) 

yielding W(x). The function W{x) is plotted in Figure [6j Note that the function is multivalued 
in the interval [— 1/e, 0], the upper branch (lighter shade in Figure [6]) is conventionally taken to 
be the principle branch Wq(x). However, for the expansion that we are interested in, we need the 
lower branch W~i(x) (darker shade in Figure [6]). As x approaches zero from below x — > 0_ the 
lower branch W-\{x) can be developed in the following series, 

W_i(x) « - In (-1) - In (In {-\)) - ^/"^ + • • • , (93) 
x x m[—-) 

this is the series that we will need to complete our derivation of a lower bound on the volume 
spectrum. 

Using these observations about the Lambert W function we have from (91), 

n « f ( J - J °) (94) 

W 54W_ 1 ((J-J )/(6ef 1 /2))- ^> 
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W(x) 



-i — i — i — i — r- 



_i i i i I i i i i I i i i i I i i i 

12 3 

Figure 6: The real values of the Lambert W function 



And expanding this solution for small (/ — Io)/(6ef 1//2 ) yields, 




6 ^)+ln(ln(^))+... 



(95) 



in the case that the phase space is odd dimensional, and a very similar result in the even dimensional 
case. Both cases are plotted in Figure [7] the odd dimensional case being the middle curve and the 
even dimensional case the lowest curve. Note that the smallest Bohr-Sommerfeld eigenvalues in 
this plot are in much poorer agreement than for the case of Figure [7j This is due to the fact, 



discussed in section 6.1, that the space of shapes is no longer a differentiable manifold when we 
consider equal spins ji, i.e. equal areas A[. 



7 Conclusion 

At the Planck scale, a quantum behavior of the geometry of space is expected. Loop gravity provides 
a specific realization of this expectation: it predicts a granularity of space with each grain having a 
quantum behavior. In this paper we have presented a new independent road to the granularity of 
space and the computation of the spectrum of the volume. The derivation is based on semiclassical 
arguments applied to the simplest model for a grain of space, a Euclidean tetrahedron, and is 
closely related to Regge's discretization of gravity and to more recent ideas about general relativity 
and quantum geometry |36j. The spectrum has been computed by applying Bohr-Sommerfeld 
quantization to the volume of a tetrahedron seen as an observable on the phase space of shapes. 
We briefly summarize our results: 



i. Spectrum. We studied volume orbits on the phase space of a tetrahedron. The main result 
is a closed formula for the action integral I(q) in terms of elliptic functions, eq. (67). This 
formula is then used to compute the Bohr-Sommerfeld levels of the volume. In Figures |3| 
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Figure 7: This figure compares the semiclassical scaling of the largest and smallest volume eigen- 
values (dark lines) to the Bohr-Sommerfeld (dots) and loop gravity (circles) spectra. The four spins 
are equal, j\ = 22 = j'3 = 34 = jo , and the corresponding phase space contains flat configurations. 
This explains the poorer agreement of the Bohr-Sommerfeld spectra at small eigenvalues. 



[4] and [7] and in the tables in Appendix [C] we compare the volume levels to the spectrum 
computed in loop gravity and find a remarkable quantitative agreement. 

ii. Degeneracy. Non-vanishing eigenvalues of the volume are twice generate. From the semiclas- 
sical perspective, this is understood as a consequence of the fact that there are two closed 
orbits with the same volume. Tetrahedra on the two orbits are related by parity. Quantized 
orbits with zero volume appear only when the symplectic area of phase space is an odd mul- 
tiple of the elementary, phase-space Planck cell. In this case parity sends the orbit into itself 
and the associated eigenvalue is non-degenerate. 

In section [5] we also identified a new degeneracy that connects multiplets in different in- 
tertwiner spaces. At the semiclassical level this symmetry originates from scissor-congruent 
auxiliary tetrahedra, and at the quantum level from the Regge symmetry of the {6j} symbol. 
We conjecture that the associated geometric tetrahedra are scissor-congruent as well. 

iii. Maximum volume. At the classical level, fixing the four areas Ai, . . . , A4, there is a maximum 
volume the tetrahedron can attain. In section [6T2] we computed the largest eigenvalue of the 
volume by expanding the action integral I{q) around the maximum classical volume. The 
largest eigenvalue is smaller than the maximum classically-allowed volume and the difference 
between the two is given by n divided by the period T of the volume orbit. Moreover the large 
eigenvalues are equispaced with a separation 2tt/T. This phenomenon can be understood as 
a manifestation of Bohr's correspondence principle. 

If instead of fixing the area of each of the four faces we fix only the area of three faces, then the 
tetrahedron of maximum volume is tri-rectangular and the maximum eigenvalue is given in 



equation ( 83 ) . We also considered fixing only the total surface area. In this case the maximum 
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volume is attained by a tetrahedron with faces having all the same area and the maximum 
eigenvalue is given in equation (82). This formula reproduces the scaling u max ~ j 3 ^ 2 found 



by Brunnemann and Thiemann, and determines corrections to it. 

iv. Minimum volume. We studied the volume gap for a quantum tetrahedron by pushing the 
Bohr-Sommerfeld approximation into the deep quantum regime. The elliptic parameter m <C 



1 that controls the accuracy of the approximation is given by equation ( 88 ) . For phase spaces 
containing no flat configurations, we find a volume gap v m i n = c {AiA^A^A^) 1 ^^ with c equal 
to, 2/3 for an odd and v2/3 for an even, number of levels. 

We have also studied phase spaces containing flat configurations. An example is the equi- 
area case. This situation is more delicate as there are orbits of infinite period corresponding 
to singular points in phase space. Physically this is also the most interesting situation as 
it leads to the smallest attainable volume when the areas are all equal to the smallest area 
A = 1/2. All our results support the existence of a volume gap for the 4-valent case given in 
the equi-area case by the unusual Log series of eq. (95), see Fig. [7J 



The remarkable quantitative agreement of the spectrum calculated here and the spectrum of 
the volume in loop gravity lends further credibility to the structure of this theory. The semiclassical 
methods of this paper provide a new understanding of many aspects of the rich structure of the 
volume spectrum in loop gravity and the explicit formulas open new avenues for analytical explo- 
ration. This is important because a deep understanding of the spectra of geometrical operators 
provides fertile ground for developing phenomenological tests of loop gravity |37j . 
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A Derivation of volume matrix elements 

The simplest derivation of the matrix elements of the volume operator that we know of is also the 
oldest and is due to Levy-Leblond and Levy-Nahas, [27]. The presentation in this Appendix closely 
parallels their argument. 
Note first that 

[Ji • J 2 , Ji ■ J 3 ] = [Ju, Jij}44 = ieij k Ji44 = Mi ■ x -h), (96) 

which is also \J\ ■ J%, J\ • J 3 ] = \ [{J\ + J2) 2 , Ji • Jz]- This allows one to express the matrix elements 
of Q by 

(k\Q\k') = Q k y = (fc|Ji • (J 2 x J 3 )\k') = -Uklhli + J 2 ?Ji • f 3 ]\k>) 

2 2 (97) 
= - l -(k(k + 1) - k'(k! + l))(k\Jx ■ J 3 \k'), 



where the last equality follows from our definition of the basis (Eq. (13) and above). Already 
it is clear that the diagonal matrix elements vanish. Then, the problem has been reduced to 
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evaluating the matrix element (k\J± ■ Js\k') for k ^ k! . Evaluating this matrix element amounts 
to two applications of the Wigner-Eckart theorem, which are performed here graphically. The 
graphical notation is quite efficient as long as you do not track phases. The overall phase is not 
needed anywhere in this work and so we proceed without tracking phases, our end result is in 
agreement with [27]. 

Inserting two resolutions of the identity yields, 



\Jl • J?>\k') = ao(k,k) x 



A 



3i 



Ji 



V 



k' 

h h 



1 

k J3 



V 



(98) 



where the constant ao is given by, 

a (k, k') = y/W+ly/W+ly/hih + Wh + VVW* + !)( 2 i3 + !)• (99) 
Graphs separable on three lines can be simplified, [38J, and in this case the reduction yields, 



{k\J 1 -J 3 \k')=a Q (k,k') x 





(100) 



Each of the simplified graphs is a 6j-symbol and it immediately follows that the matrix elements of 
the volume operator satisfy selection rules: Each of the nodes at the center of the diagram imposes 
the condition that the matrix element vanishes unless k and k 1 differ by one. This implies that the 
matrix is of the form 

/ iai ••• \ 



-ia\ 




ia>2 
-ia,2 



ion 



here the real are defined by, 



a k -i=iQ k k ~ l (k = 2,...,d+l), (102) 
where, once again, d is the dimension of the intertwiner space. Combining the results from (|97|) 
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and (100) it follows that the are given by, 
1 



-(k(k + 1) - (jfe - l)jfe)ao(fc, fc - 1) 



fc- 1 1 fe 



1 fe — 1 k 

ji fa fa 



^ 1 V(h +32 + k + l)(ji + fa - k + l)(ji - j 2 + fc)(j2 ~ Ji + g (103) 
4 V2/C + 1 

x \/(j3 + J4 + fc + l)(j 3 + j 4 - k + l)(j 3 - J 4 + fc)(j 4 ~ fa +g 

X - 1 

B Root Power Series 

The roots r% (i = 1, . . . , 4) are defined as the solutions to the quartic (in A 2 ) polynomial equation, 
P(A 2 ,Q 2 ) = (4A) 2 (4A) 2 - (2.4) 2 (9Q) 2 = P {A 2 ) - (2A) 2 (9Q) 2 = 0, (104) 

where, 

P (A 2 ) = [A 2 - (A 1 - A 2 ) 2 ][A 2 - (A 3 - A^A, + A 2 ) 2 - A 2 ][(A 3 + A 4 ) 2 - A 2 ]. (105) 

The barred roots fj (i = 1, . . . , 4) are defined as solutions to Pq{A 2 ) = 0. 
For the case where the fi are distinct we have, 

n = fi + hQ 2 + &AfQ 4 + . . . , (106) 

where 

K = n _ s and & = I I-V^-^- | , (i,j = l,...,4). (107) 

When fi and f2 coincide, say at f, the above series are singular and are instead replaced by, 

n = f±[iQ + vQ 2 ± P Q 3 + aQ i ±--- (i = 1,2), (108) 
the lower signs for r\ and the upper for r% and with, 

,= - i8 ^ , ^j i8)2( r- r "\, dog) 

V / (^3-r)(f 4 -f)' 2(f 3 -f) 2 (f 4 -f)2' ^ ' 

_ (18) 3 (r§r| + 4rr 3 r 4 (r 3 + r 4 ) - 14r 2 r 3 r 4 + 5f 4 ) _ 
8\/f (^3 — f) 7 / 2 (f 4 — f) 7 / 2 
(18) (f|f 2 (f3 + f 4 ) + ff3f 4 (f3 — f 4 ) 2 — 5f 2 f3f 4 (f3 + f 4 ) + 10f 3 f3f 4 — 2f 5 ) (111) 

2(f 3 — f) 5 (f 4 — f) 5 

If f 3 and f 4 coincide, say at f, the above formula holds with the replacements F3 — >• Fi, F 4 — )■ f2 
and f —¥ f. There is one more set of singular cases: when fi is zero then r% is zero for all Q and 
the series expansions of the other roots changes. For F2, F3 and f 4 distinct these are, 
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This series is well behaved when r 2 is zero, i.e. when n = r 2 , but singular when r% = = r and 
so we have one final case when n = and the larger roots coalesce: 

r 2 = f 2 + 7 i%Q 2 + 7 ^%Q 4 + ..., (113) 
(r — V2) (r — r 2 J° 

n = r k ± -fS=Q - ^%- 2 Q 2 ± J (1 % /2 Q 3 - r^Q 4 • • • (* = 3,4), (H4) 

Vr-r 2 2(r - r 2 ) 2 8(r-r 2 )'/ i (r - r 2 ) t ' 

with the lower signs for and the upper signs for r&. 



C Tables 

The tables below compare the Bohr-Sommerfeld and loop gravity spectra systematically, beginning 
with the smallest allowed spins and considering all allowed spins up to (ji, 02,33,34) = {\, \ 3). 
To illustrate the rapid improvement of the Bohr-Sommerfeld approximation, as well as the Regge 
symmetries, two higher spin cases are also included. 



Table 


31323334 


Loop gravity 


Bohr- 
Sommerfeld 


Accuracy 


1111 
2222 


0.310 


0.252 


19% 


1113 
2 2 2 2 








exact 


2 2 1 1 


0.396 


0.344 


13% 


HlO 
2 2 1 L 








exact 


113 3 
2 2 2 2 


0.463 


0.406 


12% 


1 1 1 3 

2 1 1 2 


0.498 


0.458 


8% 


1111 




0.620 




0.565 


exact 
9% 


113 5 
2 2 2 2 








exact 


1 1 1 5 

2 1 1 2 








exact 


2 2^2 


0.521 


0.458 


12% 


1 1 3 

2 1 2 Z 


0.577 


0.535 


7% 


1112 


0.620 


0.597 


4% 


13 3 3 
2 2 2 2 


0.620 


0.597 


4% 


1 1 3 3 
1 1 2 2 




0.752 




0.706 


exact 
6% 
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7i 7o 7q 1a 


Loop gravity 


Bohr- 
Sommerfeld 


Accuracy 


Ho'? 

2 2 Z 








exact 


1 1 3 o 

2 1 2 








exact 


1113 








exact 


115 5 
2 2 2 2 


0.573 


0.504 


12% 


1 195 

2 1 Z 2 


0.644 


0.599 


7% 


1 *3 Q c; 

1 O O 

2 2 2 2 


0.664 


0.621 


6% 


11 3 5 
1 1 2 2 


0.713 


0.692 


3% 


- - 2 2 

2 2 Z Z 


0.713 


0.692 


3% 


112 2 








exact 


0.858 


0.812 


5% 


1 3 3 n 
1 2 2 Z 




0.903 




0.867 


exact 
4% 


3 3 3 3 

O O O Q 
Z Z Z Z 


0.537 
0.992 


0.452 
0.947 


16% 

5% 


115 7 

2 2 2 2 








exact 


2 1 Z 2 








exact 


13 3 7 
2 2 2 2 








exact 


11 3 I 
1 1 2 2 








exact 






0.620 




0.546 


exact 
12% 





1.828 


1.795 


1.8% 




3.204 


3.162 


1.3% 


6 6 6 7 


4.225 


4.190 


0.8% 


5.133 


5.105 


0.5% 




5.989 


5.967 


0.4% 




6.817 


6.799 


0.3% 




1.828 


1.795 


1.8% 




3.204 


3.162 


1.3% 


11 13 13 13 


4.225 


4.190 


0.8% 


2 2 2 2 


5.133 


5.105 
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5.989 
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